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[ij  A  dual-beam  interferometric  synthetic  aperture  radar  measures  remotely  two  radial 
components  of  the  ocean  surface  current  from  a  single  flight  pass.  Combining  two  passes 
over  the  same  area,  all  three  orthogonal  components  of  the  surface  velocity  can  be 
retrieved.  An  experiment  is  conducted  near  the  Gulf  Stream  (GS)  boundary.  A  sharp 
change  of  the  surface  velocity  of  about  1  m/s  over  a  500  m  lateral  distance  is  measured. 
The  wind  and  wave  condition  is  dominated  by  a  14-s  swell  system  and  low  w  ind  velocity. 
The  wave  variance  inside  GS  is  about  twice  the  wave  variance  outside  the  GS  in  the 
present  data  set.  The  difference  in  the  wave  variance  is  considerably  higher  than  that  can 
be  expected  from  wave-current  interaction.  An  ocean  current  system  with  strong  shears 
such  as  the  GS  is  a  wave  guide  and  can  trap  waves  with  the  right  combinations  of 
wavelengths  and  propagation  directions.  Numerical  calculations  suggest  that  the  wave 
properties  of  the  data  set  may  satisfy  the  conditions  of  wave  trapping  by  the  GS.  The 
standing  wave  pattern  on  the  GS  side  of  the  sharp  velocity  front,  indicative  of  the  long 
swell  bouncing  off"  the  current  front,  also  offers  support  for  the  wave  guide  hypothesis.  In 
this  respect,  the  Gulf  Stream  can  be  considered  the  nature's  hydraulic  breakwater  that  can 
attenuate  about  50%  of  the  incident  wave  energy  generated  by  storms.  Its  role  in 
protecting  the  U.S.  coastlines  in  the  Atlantic  Ocean  cannot  be  overstated. 
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L  Introduction 

[:]  An  a  long-track  interferometric  synthetic  aperture 
radar  (InSAR)  can  measure  the  ocean  surface  current 
remotely  by  detecting  the  phase  difference  of  the  radar 
returns  from  the  same  surface  roughness  using  two  receiv¬ 
ing  antennas  mounted  on  an  aircraft  or  satellite  {Goldstein 
and  Zehker,  i9K7;  Goldstein  et  al,  1989;  Graher  et  at., 
1996).  The  phase  difference  between  the  received  signals 
from  the  two  antennas  is  produced  by  the  Doppler  fre¬ 
quency  shift  caused  by  the  ocean  current  advecting  the  short 
surface  waves  that  scatter  the  radar  signals  back  to  the 
receivers.  To  the  first  order  of  approximation,  the  phase  shift 
is  proportional  linearly  to  the  velocity  component  projected 
in  the  range  direction.  With  proper  design,  the  InSAR  can 
yield  sufficient  spatial  resolution  for  surface  wave  measure¬ 
ment  [Murom  et  al. ,  1990.  1991:  Shemer  et  ai,  1993; 
Goldstein  et  al .,  1994}.  This  is  a  significant  improvement 
lor  remote  sensing  of  ocean  waves  using  the  SAR  technol¬ 
ogy  because  the  relationship  between  surface  waves  and  the 
InSAR  velocity  product  is  much  more  straightforward  than 
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the  modulation  transfer  function  associating  ihe  surface 
waves  and  the  SAR  signal  amplitude  [c.g.,  Alpers  and 
Hasselmann,  1978:  Hasselmann  et  a l ..  1985].  Because 
distortion  of  the  waveform  through  the  velocity  bunching 
mechanism  remains  in  the  InSAR  measurements,  reliable 
re  I  riev  a  1  o  f  t  he  wa  v  c  i  n  fo rmat  i  on  fro  m  I  nS  AR  is  si  1 1 1 1  imit ed 
b>  the  nonlinearity  of  the  wave  field  [e  g..  Itwhon  et  al* 
1999:  Bao  et  al .,  1999;  He  and  Alpers.  2003].  Despite  lids 
limitation,  the  development  of  InSAR  is  indeed  a  significant 
step  forward  for  remote  sensing  of  ocean  waves. 

[}]  In  the  side- looking  InSAR,  only  one  radial  compo¬ 
nent  in  the  radar  range  direction  is  measured.  To  obtain  the 
vector  field  of  ihe  ocean  current t  multiple  non-parallel 
passes  through  the  same  ocean  surface  area  have  to  be 
conducted.  The  concept  of  deriving  two  velocity  compo¬ 
nents  from  the  In  SAR  return  by  splitting  the  radar  beam  is 
discussed  by  Rodriguez  et  al  [  I995J.  Frasier  and  Camps 
[2001]  advance  the  idea  to  employing  two  pairs  of  antennas 
placed  with  two  different  squint  angles,  one  pair  looking 
fore  and  one  pair  all.  and  each  pair  serves  as  an  InSAR 
system.  The  images  from  the  two  separate  pairs  (fore  and 
aft)  are  then  geo -collocated  to  yield  the  velocity  vector  field 
of  the  imaged  ocean  surface  area  in  a  single  flight  pass,  The 
description  of  the  system,  called  a  dual-beam  interferometer 
(DB1),  has  been  given  by  Farquharson  et  al.  [2004]  and 
Toporkov  et  al.  [2005],  the  latter  paper  also  reports  an 
application  of  the  system  to  map  out  the  velocity  field 
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around  two  tidal  inlets  in  the  barrier  islands  west  of  Fort 
Myers.  Florida.  The  retrieved  current  Field  follows  the 
outflow  pattern  expected  from  the  geometry  of  the  barrier 
islands  and  the  inlets.  Comparisons  with  l he  tidal  current 
magnitudes  predicted  by  the  U.S.  National  Ocean  Service 
reveal  discrepancies  of  up  to  0.5  m/s  Their  analyses  suggest 
that  an  important  factor  contributing  to  the  discrepancies  is 
the  effect  of  ocean  surface  waves  to  the  overall  InSAR 
velocity  measurement.  Numerical  computations  indicate 
that  waves  of  a  Pierson-Moskowit/  spectrum  for  the  ob¬ 
served  wind  speed  (5  m/s  from  NHL)  and  the  inlet  geometry 
may  have  contributed  to  the  mean  Doppler  frequency  shifts 
the  equivalent  of  -0.24  ni/s  in  the  east-west  velocity 
component  and  0.33  m/s  in  the  north-south  component. 

I  a]  Because  surface  waves  ulfect  the  accuracy  of  the 
surface  current  derivation,  and  that  wave  in  format  ion  is 
contained  in  the  InSAR  data,  in  this  paper,  we  investigate 
the  computation  of  surface  wav  e  spectrum  using  the  DBl 
measurements.  The  radial  velocity  is  contributed  by  all 
three  orthogonal  components  of  the  surface  current.  In  ihe 
2D  solution  to  retrieve  the  surface  velocity  from  the  DBl 
output,  the  vertical  velocity  component  is  assumed  lo  be 
zero  in  the  data  processing  ] Frasier  and  Camps,  2001; 
Toporkav  et  cd.<  2005].  In  ibis  paper,  ihe  full  3D  solution 
to  obtain  all  three  surface  velocity  components  by  com¬ 
bining  two  flight  passes  is  derived.  In  combining  measure¬ 
ments  at  different  times,  the  ergodic  property  of  current 
field  is  assumed  implicitly  and  the  time  Jag  is  equivalent  to 
a  phase  lag,  thus  does  not  change  the  spectral  properties 
(Appendix  All. 

\y]  In  the  course  of  wave  spectral  analysis,  it  is  found  that 
for  this  data  set  (dominated  by  a  swell  system  of  14-s  period 
and  low'  wind  speed)  the  differences  in  the  wave  properties 
on  the  two  sides  of  the  GS  velocity  front  is  much  larger  than 
that  can  be  explained  by  the  current  modulation  of  surface 
waves.  Further  analysis  suggests  dial  a  more  likely  expla¬ 
nation  of  the  observed  strong  enhancement  of  the  waves 
inside  the  GS  is  due  to  wave  trapping.  That  is,  the  wave 
properties  of  the  daia  set  satisfy  the  conditions  for  the  GS  to 
serve  as  a  wave  guide.  This  hypothesis  is  further  supported 
by  the  conspicuous  standing  wave  pattern  (of  the  dominant 
swell  length  scales)  observed  in  the  GS  side  of  the  current 
from,  indicative  of  the  swell  bouncing  off  the  current 
boundary. 

[ft]  In  the  following,  section  2  describes  ihe  DBl  current 
mapping  experiment  near  the  GS  front  cast  of  Cape  Can¬ 
averal,  Florida.  Section  3  presents  data  processing  proce¬ 
dures  and  section  4  the  wave  spectral  analysis.  Section  5 
discusses  the  wave  guide  hypothesis  and  Ihe  resolution 
issues,  and  a  summary  is  given  in  section  6. 

2.  DBl  Measurements  Near  the  Gulf  Stream 
Boundary 

2,!,  Environmental  Conditions 

[  -  ]  I  n  M  a  rc  h  2  004 .  a  n  ex  pc  r  i  mem  w  us  ct  >  nd  uc  ted  t  o  i  est 
the  DBl  system.  The  general  location  of  the  experiment  is 
offshore  of  Cape  Canaveral,  Florida  (Figure  I ).  Two  NDBC 
buoys  (41009  and  41010)  are  nearby  and  provide  the 
pertinent  environmental  information  including  wind  veloc¬ 
ity,  air  and  w'alcr  temperatures,  and  wave  properties 
(Figure  2).  Buoy  41009  is  located  at  2K.50°N  804 7°W, 


Figure  1,  Map  of  the  experimental  site.  Circles  show  ihe 
two  NDBC  buoys  that  provide  the  in  silu  wind  and  wave 
data  used  in  this  study.  The  star  shows  ihe  location  of  DBl 
data. 

which  is  about  15  km  from  the  central  location  of  the  DBl 
data.  The  mean  local  water  depth  is  41.5  m.  Buoy  4 MHO  is 
located  at  28.95°N  78,47°Wf  about  145  km  from  the  DBl 
site,  and  the  mean  local  water  depth  is  872.6  m.  The  data  to 
be  presented  in  this  paper  are  collected  from  two  flight 
passes  on  13  March  2004.  The  starting  time  of  the  firsl  pass 
(eastbound)  is  00:30:22  UTC  (marked  with  vertical  dashed 
lines  in  Figure  2)  and  the  duration  of  the  data  segment  is 
104  s.  The  second  pass  (westbound)  occurred  about  10  min 
later  (starting  time  00:40:18,6  UTC)  with  a  duration  of 
106  s.  The  surface  area  of  the  DBl  coverage  in  each  flight 
pass  is  about  2.8  km  x  10  km.  The  location  of  the  DBl  data 
is  indicated  by  a  star  in  Figure  I.  The  local  water  depth 
is  about  1 00  m  (estimated  from  the  bathymetry  map  shown 
in  Figure  2  of  Zantopp  et  uL  [1987]).  Figure  2  shows  the 
time  series  of  wind  speed,  l/m,  wind  direction,  IV*  signif¬ 
icant  wave  height,  peak  wave  period,  /),.  air  1  cm  pc  na¬ 
ture,  Tu ,  and  water  temperature,  TV  recorded  from  ihe  two 
buoys.  As  shown  in  (he  figure,  a  couple  of  high- wind  events 
passed  through  the  area  in  ihe  week  before  data  acquisition. 
The  event  on  1 0  to  1 1  March  had  sustained  wind  speeds 
between  10  and  12  nvs  and  lasted  for  about  16  hours.  The 
weather  system  continued  moving  eastward  and  ai  the  lime 
oT  DBl  experiment  more  than  one  day  later,  ihe  significant 
wave  height  at  the  offshore  buoy  site  (41010)  is  almost  5  rn 
high  although  the  wind  speed  dropped  to  below  5  m/s  at 
both  buoy  locations.  The  wave  height  in  the  nearshore  site 
(41009)  is  about  one-half  of  the  offshore  magnitude.  The 
peak  wave  periods  reported  by  the  two  buoys  are  14.29  s 
(41010)  and  13.79  s  (41009)  at  the  lime  of  DBl  data.  The 
water  temperature  at  4 1 009  slowly  increased  during  the 


2  of  1 5 


envoi  4 


HWANG  ET  AL-:  INSAR  ANALYSIS  OF  (illLF  STREAM  WAVE  GUIDE 


COTON 


25 

u  20 

e 

15 


-  . -  4JQ09  water  * 

*■  +  ... 
\  .  ^  ^  '  V 

**  -  "  «  -  41010  water  _ 

'  -  41009  air 

. 

v  *  „ 

41010  air  - 

1  (c) 


10  11  12  13  14  15  16  17 

Time  (UTC  Mar  2004) 


Figure  2,  Time  series  showing  (a)  wind  speed,  (bt  wind  direction,  (e)  significant  wave  height.  (d)  peak 
wave  period,  and  (e)  air  and  water  temperatures  measured  by  buoys  4 10(b)  and  41010.  The  vertical  dotted 
line  in  each  panel  indicates  the  time  of  DBI  data  acquisition. 


week  prior  to  the  DBI  data  acquisition  and  reached  the  same 
temperature  as  that  at  41010  on  10  March  2004,  suggesting 
that  the  (iS  was  moving  shoreward  during  l hat  period, 
f  igure  3  shows  the  2.55-day  average  sea  surface  tempera¬ 
ture  map  provided  by  the  Ocean  Remote  Sensing  Group, 
Applied  Physics  Laboratory,  Johns  Hopkins  University 
(http;  7srbdaiajluiapl.edu/d0043/a\  hrr;gs/averages  04mar' 
gs  04mar13  2016  multi. png).  The  locations  of  the  two 
buoys  arc  sketched  on  the  map  for  reference;  they  are 
approximately  on  the  two  opposite  ends  of  the  thermal 
boundary  depicting  the  G$. 

2,2.  Instrumentation 

[*|  A  deluded  analysis  of  the  DBI  design  Is  given  by 
h' raster  ami  Camps  [2001  ].  The  DBI  configurations  used  in 
the  present  experiment  have  been  reported  by  Farquharson 
etal.  (2004)  and  Tnporkav  et  af.  [2005]  A  brief  summary  is 
given  here.  The  DBI  operates  at  C-band  (5.3  GHz)  and 
vertical  polarization,  h  emits  a  6.25-ps-long  chirp  signal 
with  a  25-MI  1/  bandwidth  (hat  provides  a  6-m  range 
resolution.  The  squint  angles  of  the  fore-  and  all- looking 
antenna  pairs  are  nominally  +20"  and  -20°,  respectively. 
The  physical  baseline  in  each  pair  is  1,23  m,  only  the  fore¬ 
looking  antenna  in  each  pair  transmits.  The  antennas  point 
at  a  70°  incident  angle  in  their  squinted  planes.  The  antenna 
patterns  are  broad  in  elevation  (31°J  and  narrow  in  azimuth 
(7°)  The  aircraft  speed  is  nominally  100  m/s  and  the 
altitude  600  m.  The  range  of  the  incident  angles  of  the 
image  pixels  is  between  50  and  SI*.  Figure  4  shows 
the  phase  diagrams  derived  from  the  fore-  and  a  ft- looking 
InSAR  pairs,  Positive  v  is  in  the  direction  of  the  ftight.  The 
DBI  is  mounted  on  the  port  side  of  the  aircraft  and  the 


projected  range  increases  in  the  positive  v  direction.  1  he  top 
panel  is  the  phase  map  derived  from  the  aft-looking  pair  of 
antennas  and  the  lower  panel  from  the  fore-looking  pair. 
Image  boundaries  are  tilted  because  of  the  antenna  squints. 
Noisy,  horizontally  striped  margins  on  left  and  right  of  each 
image  lie  outside  the  synthesized  data  range  and  contain  no 
physical  information.  Significant  change  of  the  phase  pro¬ 
duced  by  the  GS  (on  the  right-hand  half  of  the  image)  is 
clearly  visible  in  both  images.  In  the  following  discussions, 
the  overlap  region  of  the  two  images,  marked  by  two 
vertical  lines  in  each  map,  is  further  processed  to  obtain 
the  ocean  surface  current. 

[9]  Data  from  the  second  pass  are  processed  in  a  similar 
fashion.  The  common  area  (1.7  km  x  7.K  km)  of  the  two 
passes  is  then  extracted  for  deriving  the  3D  solution  of  the 
three  orthogonal  surface  velocity  components.  The  proce¬ 
dure  is  described  in  the  next  section, 

3.  Data  Processing 

[  to]  The  radial  velocity  components,  to  and  u2.  in  the  two 
squint  directions  are  related  to  the  phases,  and  ■K,  by 
[Frasier  and  Camps,  2001:  Toporkov  et  ai,  2005] 


w  here  X  is  the  radar  wavelength,  ip  the  platform  speed,  and 
Bv  the  effective  baseline,  which  is  one-half  of  the  physical 
a  long-track  antenna  separation  because  only  one  antenna  in 
each  pair  is  transmitting.  With  the  present  configuration,  the 
range  of  current  velocities  without  wraparound  ambiguity  is 
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Figure  3*  The  average  sea  surface  tempera  lure  on  13 
March  2004  processed  by  the  Remote  Sensing  Group* 
Applied  Physics  Laboratory.  Johns  Hopkins  University 
(htlp:  /srbdata jhuapi.edu/d0043/avhrr  gsavenigcs/Q4mar 
gs  04marl3_2O16  multipng).  The  locations  of  the  two 
buoys  are  sketched  on  the  map  for  reference. 


±2.30  m/s  for  the  full  range  of  the  phase  angles  (±n).  The 
raw  processed  phase  data  contain  an  arbitrary  offset,  which 
can  be  removed  if  there  are  fixed  objects  in  the  images  for 
reference.  This  is  not  the  ease  in  the  present  situation.  The 
a  ire  rail  makes  repeated  passes  over  the  region;  the  time 
interval  between  two  consecutive  passes  is  about  600  s. 
Feature  tracking  is  applied  to  estimate  the  velocities  of 
several  slick-like  features  in  the  radar  scatter  amplitude 
maps  between  two  consecutive  passes.  The  features  in  the 
inshore  side  of  the  GS  usually  maintain  their  coherent 
characteristics  during  the  two  passes  and  are  easy  to 
identify.  Features  inside  the  GS,  however,  are  distorted 
beyond  recognition  between  two  passes  and  velocity 
estimates  from  feature  tracking  arc  only  available  Tor  the 
inshore  side  of  the  GS  front.  The  average  velocity  amplitude 
and  direction  of  five  identifiable  features  are  0.K2  nvs  and 
103*  (referenced  to  cast).  This  average  velocity  is  used  in 
the  data  processing  to  determine  the  offset  of  the  measured 
radar  phase. 

[n]  Following  the  notations  and  geometry  defined  in 
Figure  5  (reproduced  from  Figure  I  of  Toporhn  et  ai 
[2005])  and  assuming  that  the  current  is  confined  to  the 
horizontal  plane,  the  surface  velocity  components  can  be 
calculated  from  the  two  radial  components 

tfj  cos 0,2  -  tt2 cost), | 

“  sin(*k,  -  8,2 )  ^ 

;/2  sin  if  |  -  rjj  sin  (1,2 
sin($*j  —  fljrijsmft 

[12]  The  2D  solution  (equation  (2»  assumes  that  the 
vertical  velocity  component  is  negligible*  that  is*  v-  -  0. 
Tliis  is  a  reasonable  assumption  lor  obtaining  large-scale 
mean  surface  current  and  averaging  over  many  pixels  to 
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Figure  4.  The  phase  maps  of  the  DB1  data,  (lop)  aft-look,  (bottom)  fore-look.  Vertical  lines  in  the  two 
maps  show  the  overlap  area  used  in  the  subsequent  analysis. 
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Figure  5.  T  he  system  of  coordinates  and  symbols  used  in 
this  paper  (f  igure  I  of  Toporkov  vt  ui  [2005]).  f  ile  aircraft 
moves  from  right  to  left  ith  velocity  Vpy  the  antennas  point 
toward  the  port  side  of  the  aircraft,  one  pair  looking  forward 
with  a  squint  angle  H, the  other  pair  backward  with  a 
squint  angle  0,2’  A  common  target  on  the  ocean  surface  with 
velocity  r  is  shown  in  the  shaded  ovals*  The  range  ts  R\  and 
R2s  respectively,  for  the  fore-  ami  all-looking  antenna  pairs. 
The  altitude  of  the  aircraft  is  H  and  the  incident  angle  0,. 


suppress  the  surface  wave  effect  is  acceptable.  For  wave 
analysis,  however,  the  assumption  of  v.  0  is  no  longer 
desirable.  To  derive  all  three  orthogonal  velocity  compo¬ 
nents  from  the  DBl  measurements,  we  use  the  full  forma  - 
I a i io n  o f  the  ra< I i a  I  vel oc t ly, 

Uf  r,  sin  4-  vy  eosO^sin  -  r  cos  cos  0^.  (3) 

where,  /  -  1  and  2  fur  the  two  beams  in  each  pass.  To 
combine  the  second  llighi  pass  over  the  same  region  to 
formulate  three  equations  for  the  three  unknowns  in 
equation  (3),  the  relative  llight  direction  needs  to  be 
considered.  Explicitly,  combining  two  passes  using  the  first 
pass  as  reference,  with  both  beams  from  the  first  pass 
f  incident  angle  H,,  Nf2)  and  one  beam  from  the  second  pass 
(flight  direction  relative  to  the  first  pass,  incident  angle 
HVh  the  equation  is 


ri| 

sin  On  cos  Hi  sin  Hi  cos  Hi  eosRi 

V 

it2 

= 

sm  ( cos  flj.2  sin  fla  -  cos  l)<:  cos  l J, 2 

V, 

Hu  b}2  bo 

V 

where  sin  B^cos  -  cos  sin  If^sin  otj,  b\2  = 
H*sin  ivi  cos  0V*  sin  0neos  and  =  cos  0vl  sin  H*. 
Equation  (4)  can  he  written  as  u  B  v.  The  solutions  of  the 
three  surface  velocity  components  are 


P(i\  .  ,  _  ,  _  /  V> 

Pit  Oft  Oft 


(5) 


where  DH  is  the  determinant  of  matrix  B.  and  DHi  the 
determinant  of  the  matrix  formed  by  replacing  the  >th 
column  of  B  with  vector  w.  Additional  discussions  of  the  3D 
solutions  of  InSAR  measurements  for  both  side  loo  king 
(5L1)  and  dual-beam  (DBl)  systems,  and  a  comparison  of 
the  2D  and  3D  solutions  of  the  SL1  and  DBl  systems  arc 
given  by  P.  A.  Hwang  ct  al.  (unpublished  manuscript, 
2006f 

[i.t]  Figure  6  shows  the  three  velocity  components  of  the 
DBl  measurements  combining  two  passes.  The  mean  flow 
of  the  G$  system  in  this  region  is  primarily  northbound. 
From  in  situ  current  measurements,  the  daily  average  of  the 
northbound  velocity  at  the  core  of  the  GS  is  about  1.9  nvs 
(estimated  from  the  contour  map  in  Figure  1 3b  of  Zantopp 
etai  |  I9R7|).  The  mean  velocity  of  vv  displayed  in  Figure  6 
(middle  panel)  represents  an  instantaneous  snapshot  of  the 
surface  eurrenl  over  a  sizable  region  (the  duration  of  the 
ov  erlap  data  is  about  78  s,  the  area  of  coverage  about  1.7  km 
x  7.8  km),  the  maximum  magnitude  is  about  2  m/s,  A 
strong  gradient  in  \\  is  clearly  shown.  The  velocity  gradient 
near  the  GS  front  is  about  1  m/s  over  a  lateral  distance 
of  500  m,  corresponding  to  a  strong  velocity  shear  of  2  x 
10  1  s~l.  The  velocity  contrast  in  v*  or  v„.  is  much  weaker 
and  the  boundary  of  the  GS  is  barely  discemable  {top  and 
bottom  panels,  Figure  6). 


4,  Wave  Spectral  Analysis 
4,1,  2D  Spectrum 

[w]  With  velocity  data  like  Figure  6,  it  is  straight  forward 
to  compute  the  wavenumber  spectrum  of  the  surface  waves 
using  2D  fast  Fourier  transform  (FFT).  \taram  ct  al  |  1990. 
1991]  compare  the  wave  properties  derived  from  InSAR 
and  in  situ  pressure  gauge  array  in  southern  California.  The 
wind  speed  is  low  (2  m/s)  and  the  wave  Held  swell- 
do  minuted.  They  show  that  the  InSAR  and  buoy  data  are 
in  v  ery  good  agreement  for  waves  propagating  in  the  range 
direction  and  the  quality  of  agreement  deteriorates  for 
waves  traveling  in  the  azimuth  direction.  Similar  conclusion 
is  reached  by  Goldstein  ct  al.  [1994],  The  deteriorating 
agreement  for  wave  components  approaching  azimuth  di¬ 
rection  (with  wavelengths  much  longer  than  that  affected  by 
velocity  bunching  mechanism)  is  a  fundamental  problem  of 
wave  measurements  by  InSAR.  The  radial  velocity  is 
contributed  by  two  velocity  components  in  the  vertical 
and  ground  range  directions*  An  azimuth  traveling  wave 
gives  no  contribution  to  either  component.  Although  in 
principal,  this  is  taken  into  account  in  the  sin4>  dependence 
of  the  transfer  function  between  the  spectra  of  surface 
velocity  and  surface  elevation,  the  function  has  a  second- 
order  singularity  behav  ior  [U\/AT  -  1/sin^o,  where  k  is  the 
surface  wavenumber  and  o  the  wave  propagation  direction]. 
The  noise  in  the  system  is  greatly  amplified  for 
approaches  zero  when  the  transfer  function  is  applied  to 
convert  the  velocity  spectrum  to  displacement  spectrum. 
This  singularity  behavior  is  masked  by  the  incident  angle  (U) 
dependence  in  the  transfer  function  between  the  radial 
velocity  and  the  surface  displacement  spectra  (sin'tt  sirr<> 
+  cojtO)  1  [eg.,  llas&elmam  and  Hasselmamu  1991]. 
When  possible,  il  is  much  more  straight  forward  using  the 
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Figure  h.  Maps  showing  the  three  orthogonal  components  of  the  surface  velocity  field,  (a)  v (b)  \\*  and 
(c)  v-*  computed  from  combining  two  passes  of  the  DBI  measurements. 


surface  velocity  components  (instead  of  the  radial  velocity  ) 
for  wave  analysis,  to  be  further  explained  in  the  following. 

[  \  f|  Of  the  three  components  of  the  wave-induced  orbital 
velocity,  v,,  v,  and  vr,  the  vertical  component  is  not 
dependent  on  the  wave  propagation  direction,  as  shown 
below.  For  a  given  spectral  wave  component  (denoted  by 
subscript  a),  the  surface  displacement  qr„  the  stream  func¬ 
tion  of  the  wave  field  and  the  orbital  velocity  v„  =-  (vv„, 
vvtt*  v,tf)  arc  [e.g.,  Phillips,  1977] 

lit  -  cos{*„,v  - 


from  l  he  sum  of  the  two  spectra,  Srt  and  5,^  computed  from 
\\  and  tv  respectively,  or  from  Sy-  computed  from  That 
is,  we  can  express  the  velocity  components  as  the  Fourier 
series 

V  *».  V 

V*  -  ^  ^  l\j|4  %\  =  £***  -  ^  ^  Vztf  (^) 

H  -  I  A  I  q  I 

The  directional  spectrum  of  a  given  wavenumber  compo¬ 
nent,  4>nK  can  be  written  as 


9m 


4 Huu  cosh  Ay[r  4-  h) 
k„  sinh  kji 


sinfAji  t 


~4«f  4  £„). 


/  rArni  *  *'rn  ^£m  ^ 
\  (hr  fh  ih  } 


m 


-  eos :o„5v(t 

•^n  (^Vr .  )  “  SIH  ■  (9) 

Svi  .  Ojij)  *V  (An „  t’y) 


vv  here  A  is  the  wavenumber  vector  with  modulus  A,  a  wave 
amplitude.  4  angular  frequency,  h  water  depth,  /  time  and  - 
a  random  phase,  l  or  InSAK  applications,  £  0  (for  velocity 

at  the  surface)  and  /  can  be  set  to  zero  without  loss  of 
generality.  Explicitly,  the  three  components  of  the  surface 
velocity  are 

v\„  Ait** p  cost  k„x  -  rh  !  eos  o„ 

vVm  -  cost kflx  4  cw  )  sin  4»if  .  (7) 

V2tt  -  sin(^MAT  4  „  1 

The  surface  displacement  spectrum  Su  can  be  calculated 
from  the  velocity  spectrum  S,.  The  latter  can  be  derived 


The  surface  displacement  spectrum  is  related  to  the  velocity 
spectrum  by 

Sn(k}  =  uT3£<*)  (JO) 

The  wave  angular  frequency  is  related  to  wavenumber  by 
the  dispersion  relation,  a2  -  gk  tanh  Ml  where  g  is  the 
gra\  ilational  acceleration. 

[jr*]  The  region  near  the  Gulf  Stream  boundary  is  dy¬ 
namically  active.  Despite  extensive  efforts  in  the  processing 
of  motion  compensation,  considerable  platform  motions 
remain  in  the  resulting  velocity  product.  This  is  detectable 
visually*  for  example,  as  the  low-frequency  undulations  m 
the  phase  and  velocity  maps  (e.g..  Figures  4  and  6f  For 
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Figure  7.  Comparison  of  Ihe  2D  wave  spectra  (a)  inside  and  (b)  outside  the  Gulf  Stream,  (c)  1 D  spectra 
from  DBI  and  buoy  measurements. 


large-scale  currents  such  as  the  GS  or  tidal  flows  through 
inlets,  these  residual  errors  of  the  platform  motion  can  be 
reduced  significantly  by  spatial  smoothing.  For  the  scale  of 
surface  gravity  waves,  the  residual  errors  in  the  platform 
motion  correction  become  a  serious  contamination  to  the 
wave  signal  As  analyzed  in  section  5,  the  signal  to  noise 
ratio  (SNR)  of  the  wave  data  (in  terms  of  the  variance  ratio 
of  waves  and  measurement  uncertainties)  is  about  two  to  six 
and  represents  a  challenge  in  wave  spectral  analysis.  An 
empirical  scale  decomposition  (BSD)  procedure  similar  to 
(lie  empirical  mode  decomposition  method  of  Huang  ei  uL 
|1998,  1999]  is  developed  to  decompose  the  measured 
signal  into  several  components,  each  with  a  narrow  band 
of  length  scales.  The  spectral  analysis  is  performed  on  each 
signal  component  to  reduce  the  contamination  of  spectral 
signal  in  the  wavelength  band  by  the  spectral  leakage  of 
lunger-scale  waves,  which  arc  usually  of  much  higher 
spectral  density  the  wavenumber  dropoff  of  a  typical 
wind  wa\e  spectrum  is  k  l\  More  detail  on  the  BSD  is 
given  in  Appendix  A2, 

[it]  Figures  7a  and  7b  display  the  2D  velocity  spectra, 
k\l  inside  and  outside  GS  calculated  with  v-.  As 
described  in  the  last  section,  the  wave  field  is  the  remnant 
from  an  earlier  northerly  event  (Figure  2).  This  information 
can  be  used  to  resolve  the  1 80°  ambiguity  in  the  2D  spectra 
and  u  is  reasonable  to  assume  that  the  dominant  directions 
of  the  wave  components  shown  in  F  igure  7  are  from  Ihe  first 
quadrant.  Significant  wave  height  and  peak  wave  period 
(//,.  Tfr )  at  the  time  of  radar  data  are  (t.15  m,  13.8  s) 
recorded  by  Buoy  41009,  and  (2.60  m,  14.3  s)  by  Buoy 
41010;  the  corresponding  peak  wavenumbers  are  2.65  x 


10  3  and  1.97  x  10  3  rad/m,  respectively.  The  spectral 
peak  wavenumbers  derived  from  the  DBI  analysis  are  in 
reasonable  agreement  with  the  buoy  measurements  and  will 
be  further  discussed  in  Ihe  next  subsection  describing  the 
properties  of  the  ID  wave  spectrum  (Figure  7c k  Poorer 
resolution  in  the  azimuth  direction  due  to  velocity  bunching 
[e.g..  A/perS  and  Rufetnu  lu  1979;  Hasselmann  et  ai,  1985] 
is  clearly  shown  in  the  contours  of  the  2D  spectra.  For  the 
present  data  set,  the  azimuth  roll-off  becomes  serious  at 
about  ks  0,25  rad'm,  which  gives  a  Nyquist  wavelength  of 
25  m  and  the  effective  resolution  in  the  azimuth  direction  is 
12.5  rn.  The  nonlinearity  of  SAR  measurements  can  be 
quantified  by  the  dimensionless  parameter  kpu>R/l‘p  |e,g., 
Bao  ei  a/.,  1999],  w  here  ur  is  the  magnitude  of  the  wave 
orbital  velocity.  For  the  present  data,  R  2500  m,  lr  ^ 
1 00  m/s,  kp  -  6.02  rad/m,  and  u,  ^  0.22  to  0.58  m/s  for  a  swell 
of  14-s  period  and  0.6  to  1.3  m  amplitude,  the  nonlinearity 
parameter  is  about  0. 1 3  to  0.29  and  the  process  of  SA R  image 
formation  in  the  present  data  set  is  reasonably  linear. 

[is]  Waves  inside  the  G$  are  distinctively  different  from 
those  outside  the  GS  in  three  respects,  (a)  The  spectral 
density  level  is  considerably  higher  inside  GS,  by  about  a 
factor  of  two  (Figure  A2).  (b)  The  directional  distribution  of 
the  wave  spectrum  outside  the  GS  is  rotated  counter¬ 
clockwise  by  about  20  to  10  degrees  in  comparison  with 
that  of  Ihe  spectrum  inside  the  GS.  And  (c)  there  is  a 
conspicuous  standing  wave  (cross-hatched)  pattern  in  the 
dominant  scale  signal  components  in  the  G$  side  of  the 
sharp  current  front  (Figure  Al). 

[w]  The  directional  difference  of  the  two  wave  fields  is 
consistent  with  that  expected  from  wave-current  interaction. 
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This  can  he  illustrated  hy  the  wavenumber  conservation 
equation  (c,g„  sect  ion  2  6  of  Phillips  [1977]) 

Ok 

—  (II) 

Or 

where  n  is  the  apparent  frequency  of  the  wave  field,  which 
is  related  to  the  intrinsic  frequency,  a),  by  ft  -  +  A  ■  V  and 

r  the  current  vector  Assuming  that  the  intrinsic  frequency 
does  not  vary  with  space,  the  following  equation  describes 
the  evolution  of  the  wavenumber  vector  propagating 
through  a  current  (kid 


Ok 

i)t 


Vikt y 


(12) 


To  apply  this  to  the  C tS  in  the  experimental  location,  the 
current  lie  Id  is  simplified  to  V  (0,  t  \  and  dli)y  0;  thus 
(12)  becomes 


<2*|  0{k2n 

Or  Ox 


f  13} 


Figure  8  shows  an  example  of  the  trajectories  of  waves  with 
a  southbound  component  going  through  a  current  field  F(.v) 
f  tiCxp[  with  </  N,  scaling  width  a*  -  80  km,  and 

maximum  velocity  F,t  2  m/s.  (Note,  even  with  </  -  8  for 
the  velocity  profile,  the  maximum  current  shear  is  7.4  x 
10  v  s  \  si  ill  much  smaller  than  the  observed  value  of  2  x 
10  s  1  from  the  smoothed  velocity  map  of  \\  shown  in 
Figure  6).  In  this  plot,  the  reference  direction  of  wave 
propagation  is  0  degree  in  the  opposite  direction  of  current* 
The  evolution  of  wave  propagation  is  consistent  with  the  2D 
wavenumber  spectra  shown  in  Figure  7* 

4,2,  ID  Spectrum 

[:n]  The  PBI  2D  wavenumber  spectra  can  be  integrated 
to  yield  ID  spectra  to  compare  with  the  buoy  measurements 
(Figure  7c),  The  buoy  spectrum  is  in  frequency  domain  and 
the  following  equation  is  used  to  convert  the  frequency 
spectrum  to  wavenumber  spectrum, 

S(k)^.m (14) 
ak 

As  discussed  in  section  2  (Figure  2),  the  waves  at  the 
experimental  site  are  dominated  by  the  swell  generated  by 
an  earlier  northerly  event*  The  wavenumbers  at  the  spectral 
peaks  at  the  two  buoy  sites  are  slightly  different,  2.65  x 
1 0  '  and  1 .97  x  10  2  rad/m,  respectively,  for  4 1 009  and 
41010  (about  160  km  apart).  The  difference  in  the  spectral 
peak  wavenumbers  of  the  two  buoy  measurements  can  be 
explained  by  the  refraction  effect  due  to  shoaling  from  deep 
to  shallower  water:  the  water  depth  is  873  m  at  buoy  41010 
and  4L5  rn  at  buoy  41009.  Another  factor  is  the  Doppler 
frequency  shift  (not  included  in  the  frequency  to  wave¬ 
number  conversion  here).  For  waves  udvccted  by  a  counter 
current  with  projected  velocity  of  I  m/s  in  the  direction  of 
wave  propagation,  the  Doppler  frequency  shift  introduces 
about  4  percent  increase  in  the  apparent  frequency  for  k 


Figure  8,  Numerical  computations  illustrating  the  wave 
trapping  effect  of  the  Gulf  Stream  wave  guide,  k0  ^  2  x 
10  *  rad. in.  The  reference  direction  of  wave  propagation  in 
the  legend  is  0  degree  in  the  opposite  direction  of  current, 

0.02  rad/m.  Using  the  energy  (lux  conservation  principle, 
Si\(k)Cf£k)  =  constant,  the  shoaling  cl  feet  (without  external 
sources  and  sinks)  can  be  quantified,  f  igure  9  shows  the 
buoy  and  DB1  spectra  adjusted  for  shoaling  refraction  to  the 
water  depth  of  buoy  41009,  The  spectra  computed  from 
DBI  are  from  two  locations  about  2,7  km  apart  (between  the 
centers  of  the  two  squares  used  lor  spectral  analysis)  on  the 
two  sides  of  the  sharp  velocity  front,  1 5  km  from  the  inshore 
buoy  (41009)  and  145  km  from  the  offshore  buoy  (41010), 
The  re frael ion -adjusted  spectral  peak  wavenumbers  of  DBI 
and  buoy  are  ufl  near  2.7  x  10  2  rad/m. 

[21]  The  most  interesting  part  of  this  comparison  is  the 
difference  of  the  two  DBI  spectra  on  (he  two  sides  of  the 
velocity  front.  To  the  GS  side,  the  DBI  spectrum  is 
essentially  the  same  as  that  of  buoy  41010.  especially  when 
the  current  modulation  is  also  accounted  for.  Numerical 
computation  of  wave-current  interaction  shows  a  level  of 
about  10  percent  increase  in  wavenumber  and  wave  vari¬ 
ance  inside  the  GS  for  wave  components  with  initial 
(unperturbed)  wavenumber  of  0.02  to  0,03  rad/m 
(Appendix  B),  The  wave  spectrum  of  buoy  41010  adjusted 
for  both  shoaling  refraction  and  wave-current  modulation  is 
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Figure  9.  ][)  spectra  From  DBl  and  buoy  measurements  adjusted  for  the  effect  of  shoaling  refraction 

using  Lhe  depth  of  buoy  41009  as  reference.  The  current  modulation  effect  of  the  deep  water  (buoy 
41010)  spectrum  is  also  shown. 


sketched  in  Figure  9.  ]|  is  fair  to  say  that  the  difference 
between  the  DBl  spectrum  in  the  CIS  side  of  the  velocity 
front  and  the  spectrum  at  buoy  41010  is  largely  attributable 
to  refraction  and  current  modulation.  This  is  not  unexpected 
because  the  wave  condition  is  dominated  by  the  14-s  swell 
The  wave  steepness  of  the  swell  is  small,  so  dissipation  by 
breaking  or  viscosity  and  generation  by  wind  (5  m/s)  do  not 
modify  this  swell  system  over  the  145-km  distance  between 
the  two  measurement  sites.  However,  just  a  short  distance 
away  to  the  inshore  side  of  the  velocity  front,  the  wave 
spectral  density  dropped  drastically  and  a  demarcation  line 
can  be  drawn  at  the  velocity  front  delineating  two  regions  of 
very  different  surface  wave  properties. 

5.  Discussions 
5.1.  Gulf  Stream  Wave  Guide 
[22 1  Increased  surface  roughness  at  the  current  boundary 
is  a  frequently  observed  phenomenon.  Modification  of 
winds  and  waves  hy  ocean  current  systems  can  also  he 
observed  from  spacebome  altimeter  measurements.  For 
example,  Hwang  [2005]  presented  art  analysis  of  the  spatial 
variation  of  wind  speeds  and  significant  wave  heights  in  the 
Yellow  and  East  China  Seas  using  six  years  data  from  the 
TOPEX  Poseidon  (TP)  output.  The  region  has  a  distinctive 
monsoon  pattern  and  the  Kuroshio  trajectory  is  stable  due  to 
bathymetric  confinement.  In  winter  months  when  w  aves  are 
propagating  obliquely  against  the  Kuroshio,  the  spatial 
distribution  of  wave  heights  shows  a  dear  modulation  by 
the  Kuroshio.  with  a  peak  enhancement  of  25  to  50  percent 
<d'the  wave  variance  observed.  Three  factors  contributing  to 
lh  c  w  ave  e  n  banc  c  men  l  a  re  c  i  ted :  i  nc  rea  sed  { e  I  Tee  live)  vv  inti 


speed  in  the  Kuroshio  domain  leading  to  an  increase  in 
wave  generation,  hydrodynamic  modulation  due  to  wave- 
current  interaction*  and  air-sea  instability  in  the  Kuroshio 
leading  to  an  enhanced  wave  generation.  The  first  factor  can 
be  calculated  from  the  wind  speed  data  available  from  the 
I  P  output.  The  second  factor  can  be  quantified  by  compu¬ 
tation  using  the  action  density  conservation  equation.  The 
result  shows  that  each  of  the  three  factors  contributes  about 
one  third  to  the  total  enhancement  observed  in  the  TP  wave 
data. 

[2i]  The  factor-of-two  increase  of  the  wave  variance  (of 
the  long  swell)  inside  the  GS  in  ihe  present  data  set  (Figure 
A2)  is  considerably  higher  than  that  expected  from  hydro- 
dynamic  modulation,  calculated  to  be  about  10  percent 
(unperturbed  wavenumber.  ki}  -  2  x  10  ■  rad/rn)  and  12 
percent  (Al}  3  x  10  ~  rad/m)  for  the  GS.  The  computation 
of  hydrodynamic  modulation  using  the  energy  or  action 
density  conservation  equation  has  been  discussed  in  great 
detail  elsewhere  [e,g,,  Longuet- Higgins  and  Stewart.  I960; 
Keller  and  WrighL  1975;  Hughes,  1978;  Phillips,  1984; 
Thompson  and  Gasparovic,  1986;  Hwang  and  Shemdin , 
1990;  Hwang*  1999,  2005]  and  a  brief  summary  is  given  in 
Appendix  B.  As  stated  earlier,  the  waves  are  dominantly 
swell  (wave  period  14  s),  thus  effects  of  air-sea  stability 
condition  and  wave-generation  by  local  w  ind  (of  5  m/s)  can 
be  ignored.  A  more  likely  mechanism  contributing  to  the 
observed  intensification  of  the  wave  field  inside  the  GS  is 
the  “wave  guide”  effect  produced  by  the  strong  velocity 
gradient  at  the  GS  boundary  .  With  the  right  combinations  of 
wave  and  current  properties,  (surface,  internal  or  acoustic) 
waves  can  be  trapped  inside  the  GS,  This  effect  occurs 
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when  the  change  of  A i  following  equation  (13)  is  so  large 
that  Aj  changes  sign  (reversing  the  propagation  direction). 
The  formation  of  a  wave  guttle  by  a  shear  current  is  a 
delicate  combination  of  several  parameters  including  the 
velocity  shear,  the  current  width*  the  wavelength  and  the 
wave  propagation  direction  (so  the  wave  guide  formed  by 
an  ocean  current  is  probably  quite  leaky),  As  shown  in 
Figure  8  for  Ay>  -2x10  '  rad/m.  wave  components  inside 
the  current  that  propagate  at  angles  less  than  ±20°  against 
the  currenl  are  trapped.  The  spectral  peak  wavenumber 
of  the  wave  field  at  the  time  of  data  collection  is  dose  to 
2  x  10  2  rad 'm  and  the  source  of  swell  is  from  northerly 
and  northeasterly  directions,  the  wave  guide  e  fleet  may  be 
important  in  producing  the  large  enhancement  of  waves 
inside  the  GS  in  the  present  data  set, 

[: i]  Another  observation  that  provides  further  support  of 
the  wave  trapping  hypothesis  is  the  standing  wave  (cross- 
hatched)  pattern  in  the  GS  side  of  the  velocity  front.  This 
feature  is  enhanced  considerably  through  the  LSD  analysis 
(Appendix  A).  The  cross-hatched  pattern  of  the  waveform 
in  the  decomposed  components  with  length  scales  near  the 
spectral  peak  wavelength  on  the  GS  side  of  the  velocity 
front  is  suggestive  of  waves  bouncing  off  the  current  front 
(see  the  velocity  maps  of  fh  and  Jl\  m  the  let!  panels  of 
Figure  At).  Based  on  the  analysis  presented  in  this  paper, 
the  Gulf  Stream  can  be  viewed  as  the  Nature's  hydraulic 
breakwater  that  protects  the  U.S.  eastern  seaboard.  Without 
it.  the  severe  waves  reaching  the  coast  from  northeasterly 
storms  may  have  been  more  energetic  by  a  factor  of  two. 

5,2.  Measurement  Uncertainty 

[:<]  As  mentioned  earlier,  there  are  considerable  residual 
noises  caused  by  the  platform  motion.  Here  the  resolution 
issue  is  further  investigated.  Frasier  and  Camps  [2001  j  and 
Toporkov  et  at.  [2005]  present  an  error  analysis  to  quantify 
l he  accuracy  of  Dill  velocity  measurements  based  on  the 
Cramer-Rao  bound  of  the  phase  error  variance  [Rodriguez 
amt  Martin.  1992:  Rosen  et  at ..  2000]. 


where  A/  is  the  number  of  looks  and  y  ihe  coherence 
parameter  between  the  radar  returns  from  the  pair  of 
antennas  forming  the  InSAR.  Hie  variance  of  the  error  in 
the  velocity  measurements  can  be  derived  from  (1)  and  (21 
assuming  that  the  phase  noise  in  the  fore-  and  aft-looking 
inierfexograms  is  uncorrelated  ]  Frasier  a  tut  Camps,  2001; 
Toporkov  et  a/.,  2005] 

:  (  V*ii  cos2 0*2  +-  fr~:  cos2 0*i 

11  "  VW  ‘  sm*  (0t|  -(),:) 

2  /  XL.  \  *o;H  sin2 Of-2  4  sin2  i\\ 

\4itZU/  sin2  ( II*  i  -  0*2 )  sin2 

With  spatial  smoothing  of  20  x  20  pixels  for  the  surface 
mean  current,  the  accuracy  of  the  velocity  measurement  is 
ox,  -  0.032  m/s  and  otl .  -  0,013  m/s  in  the  present  data  set. 
comparable  to  that  reported  by  Toporkov  et  a!  [2005]. 
Without  smoothing.  arA  -  0.33  m/$  and  o„.  -  0.16  m/s* 
which  are  very  coarse  for  surface  wave  processing.  For 


wave  spectral  analysis,  a  3  x  3 -pixel  smoothing  (2D 
running  average)  is  applied*  resulting  m  aM  0.16  m  s  and 
aiy  =  03)69  m/s.  The  root  mean  square  amplitude  of  the 
wave-induced  orbital  velocity  computed  from  the  buoy- 
measured  wave  spectrum  is  0.24  m/s  at  buoy  41009  and 
0.41  m/s  at  buoy  41010.  The  signal  to  noise  ratio  (SNR)  in 
terms  of  the  velocity  variance  is  1.9  for  waves  near  buoy 
41009  and  5.5  near  buoy  41010.  While  these  figures  arc  not 
detrimental*  special  core  is  needed  in  spectral  analysis  to 
extract  the  wave  signal  from  Ihe  noisy  background 
(Appendix  A). 

6.  Summary 

[2ft]  A  DB1  employs  two  sets  of  InSAR  to  provide 
mapping  of  both  horizontal  velocity  vectors  of  the  ocean 
surface  in  a  single  flight  pass.  In  an  earlier  study,  it  was 
shown  that  the  resolved  velocity  fields  over  barrier  islands 
follow  the  expected  outflow  pattern  of  tidal  flows  |  Toporkov 
et  at..  2005].  In  this  paper,  DBt  data  acquired  in  the  vicinity 
of  the  GS  boundary  are  analyzed.  The  3D  solution  is 
derived  to  extract  all  three  orthogonal  surface  current 
components  from  combining  two  flight  passes  over  the 
same  region.  The  retrieved  mean  current  field  (Figure  6) 
is  in  reasonable  agreement  with  in  situ  measurements 
(section  3). 

[27]  Section  4  describes  the  wave  spectra  computed  from 
the  DBt  data.  The  contamination  by  the  platform  motion  \> 
not  fully  removed  in  the  present  dataset  due  to  ihe  bumpy 
ride  caused  by  the  unstable  stratification  near  the  GS 
boundary.  The  SNR  of  the  wave  data  is  estimated  to  be 
about  two  to  six  and  considerable  care  is  needed  for  the 
spectral  analysis,  A  new  signal  processing  technique  (BSD) 
is  developed  to  improve  the  quality  of  wave  spectral 
analysis  (Appendix  A).  The  result  is  in  good  agreement 
with  in  situ  buoy  output  The  difference  in  the  propagation 
directions  of  waves  inside  and  outside  the  GS  detected  from 
the  2D  wave  spectra  is  consistent  with  the  refraction  effect 
expected  from  wave-current  interaction.  The  wave  variance 
inside  GS  is  about  twice  that  outside  CiS  in  the  DB1  data 
although  the  spatial  separation  is  only  about  2.7  km.  I  he 
level  of  enhancement  inside  CiS  is  considerably  higher  than 
that  can  be  explained  by  hydrodynamic  modulation.  Gen¬ 
eration  of  waves  by  focal  wind  is  not  an  important  factor 
judging  from  the  lung  wave  period  (14  s)  and  low  wind 
speed  (5  m/s);  so  is  the  air- sea  stability  c fleet  and  both 
factors  can  be  ignored  in  this  case,  (l-vcn  when  all  three 
factors  are  important,  earlier  investigation  indicates  25  to  50 
percent  wave  variance  enhancement  by  the  Kuroshio  using 
the  TP  data  [I twang,  2005],)  A  more  likely  explanation  of 
the  excessive  wave  enhancement  in  the  present  data  set  is 
that  the  GS  can  act  as  a  wave  guide  and  trap  waves  inside 
when  the  right  combinations  of  wavelength,  propagation 
direction*  currenl  width  and  shear  level  exist.  Numerical 
compulations  suggest  that  the  surface  waves  at  ihe  lime  of 
data  acquisition  may  satisfy  the  selective  conditions.  The 
hypothesis  is  further  supported  by  the  appearance  of  stand¬ 
ing  wave  (cross-hatched)  pattern  on  the  GS  side  of  the 
velocity  front,  indicative  of  swell  bouncing  off  the  sharp 
velocity  front  as  a  result  of  wave  trapping.  The  wave 
analysis  of  DBI  measurements  shows  that  live  Gulf  Stream 
is  Natures  hydraulic  breakwater  protecting  the  U.S.  eastern 
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Figure  AL  (a)  An  example  showing  ihe  tSD  process.  The  velocity  maps  shown  are  the  original  data 
(fix)),  live  longer  scale  components  (fi\(x),  fu(x),  .  fi_4x))*  and  the  short  scale  component 

Notice  the  standing  wave  (cross-hatched)  pattern  in  4 ;  and  f)  ^  on  the  OS  side  of  the  velocity  front,  (b)  The 
corresponding  column- hv -column  variance  of  the  original  signal  and  decomposed  components  show  n  to 
the  left.  The  sum  of  the  variance  of  scale  components  fn  t  of}*  (dominated  by  the  swell  signals)  is  added 
as  one  curv  e  to  simplify  the  figure.  Contamination  from  residual  platform  motion  is  quite  severe  for  all 
components  in  the  neighborhood  of  azimuth  distance  -^IGOO  m. 


s  ea  board .  1 1  red  uccs  by  half  the  in  ten  si  tv  o  f  se  \  e  re  wa  v  es 
reaching  the  coast  from  northeasterly  storms. 

Appendix  A:  Additional  Notes  on  Data  Analysis 

AL  Combining  Flight  Tracks  With  Time  Lags 

| :h]  When  combining  InSAR  measurements  from  differ¬ 
ent  llight  passes,  the  current  fields  on  ihe  ocean  surface  are 
obviously  taken  from  different  times.  Assuming  that  the 
spectra  of  wave  and  current  fields  are  quasi -stationary 
(ergodicX  the  orbital  velocity  field  of  a  give  spectral 
component  for  the  Mb  pass  can  be  written  as 

i\  —  j^iCosfjG  hjK'osu 

\\  =  xacasik.x  fy)  sino  .  I.A1 ) 

v  tt  uxt$\n(kx  —  6^) 

where  i\  is  the  phase  lag  U7,X  In  combining  two  passes,  say 
F)  0.  F;  we  can  introduce  a  new  time  reference  with  an 


offset  Ar  ±h/2u\  and  the 
can  he  written  as 

velocity  field  in  the  two 

passes 

v,  =uwcos(A\i 

-  jjt* 

2 )  cos  1 1 

i;  uixrcosfit.v 

-yf  - 

b/2 )  sin  *  j  , 

(A2) 

iv  ouV/  sm(  A’v 

-  u St*  - 

b/2) 

that  is,  the  time  lag  introduces  a  phase  lag  in  the  combined 
wave- induced  flow  field  and  does  not  change  the  spectrum 
This  analysis  can  be  extended  to  combinations  with  more 
than  two  llight  passes, 

[2*]  The  above  discussion  is  applied  lo  a  Fourier  spectral 
component  of  the  surface  current  field.  For  different  Fourier 
components,  the  phase  lags  arc  different  but  the  spectrum  is 
not  alTected  because  of  the  assumption  of  random  phase 
lags  inherent  in  the  spectral  analysis  In  other  word,  the 
resulting  current  field  is  an  equivalent  realization  produced 
by  the  expected  current  spectrum.  The  radial  velocity 
measured  by  the  radar  at  a  given  location  on  the  water 
surface  is  a  linear  combination  of  the  three  orthogonal 
components  of  the  surface  velocity  at  that  location  (equation 
(3)),  therefore,  the  field  of  radial  velocity  from  combining 
different  flight  tracks  is  also  an  equivalent  realization 
produced  by  the  expected  current  spectrum,  In  the  present 
data  set,  it  is  a  good  assumption  that  the  swell  condition 
remains  ergodic  during  the  10-minute  period  between  two 
llight  passes. 

A 2,  Empirical  Scale  Decomposition 

[mi]  One  of  the  difficult  problems  in  airborne  data  pro¬ 
cessing  is  the  low  SNR  because  the  complete  removal  of 
platform  motion  is  not  easy.  This  problem  also  exists  in 
signals  that  contain  many  scales,  for  example,  short  wind- 
generated  waves  in  the  radar  Bragg  resonance  scale  embed¬ 
ded  in  longer  wind  waves  and  swell  In  addition,  due  lo 
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Figure  A2.  (a)  Comparison  of  ihc  surface  wave  spectrum  computed  from  v.,  denoted  as  S\c  in  the 
legend,  and  the  sum  of  the  spectra  computed  from  vt  and  \\m  denoted  as  SV1  +  5V1  in  the  legend,  (b)  The 
ratio  between  the  wave  spectra  outside  and  inside  the  GS. 


hydrodynamic  modulation,  the  distribution  of  short  waves 
on  longer  waves  or  in  areas  with  current  shears  is  no  longer 
homogeneous  or  isotropic,  which  is  an  important  assump¬ 
tion  in  Fourier  spectral  analysis.  As  a  result  of  such  'Tiding 
wave"  problem,  the  spectral  signature  of  short  waves  is 
frequently  overwhelmed  by  the  spectral  leakage  from  longer 
scale  wave  components  because  the  spectral  density  distri¬ 
bution  of  wind  waves  scales  roughly  with  k  '  \  Huang  ct 
at.  [1998,  1999]  pioneered  an  empirical  mode  decomposi¬ 
tion  (EMD)  designed  to  reposition  the  riding  waves  at  the 
mean  water  level  The  main  idea  is  to  find  the  trend  that  can 
represent  the  mean  local  average  so  that  riding  waves  can  be 
identified.  The  EMI)  method  uses  the  point -by -point  aver¬ 
age  of  ihc  signal  envelopes  for  the  local  mean.  The 
difference  between  the  original  signal  and  the  local  mean 
represents  a  mode  of  the  signal  The  local  mean  may  also 
contain  riding  waves,  and  the  mode  decomposition  process 
continues  until  no  riding  waves  exist  in  the  local  mean 
component.  The  associated  spectral  analysis  (Hilbert-Huang 
Transformation  HHT)  provides  superior  spatial  (tempo¬ 
ral)  and  wavenumber  (frequency)  resolution  lor  nonstation¬ 
ary  (or  inhomogeneous  in  space)  and  nonlinear  data.  The 
IIHT  spectrum  also  results  in  a  considerably  different 
interpretation  of  nonlinearity  (frequency  modulation  instead 
o f  h a rmo i \ tc  gene ra t i on )  \Hw cmg  el  a /. ,  2003],  In  the  EM D 
method,  the  length  scales  of  each  mode  are  usually  mixed 
but  the  bandwidth  of  each  mode  is  generally  much  narrower 
than  that  of  the  original  signal  The  EMD  can  be  extended 
to  2D  processing  by  decomposing  the  data  in  ID  line-by¬ 
line  (Huang,  personal  communication). 


\j\]  Here  2D  digital  filter  (filter2  function  in  MAT  LAB) 
is  used  for  signal  decomposition  of  the  3D  data  such  as  the 
velocity  maps  in  this  paper.  Because  digital  filter  specifics 
length  scales  for  high-  and  low-pass  operation,  the  process 
is  called  empirical  scale  decomposition  (LSD).  Basically,  a 
signal  /(.v)  is  decomposed  sequentially  into  a  low-pass  and  a 
h  ig  h  -  pass  go  m  po  n  en  t , 

/'( vj  -/n(.v)  +  fm(j} 

Jfoifa)  =Jh  fa)  +\fc(A) 

fm  fa)  =  fu(x)  +.///*(*)  <  A3 ) 


fits  - 1  fa)  -  its  fa'  l  ffn\  faT 

which  can  also  be  written  as 

/‘fa)  =  Alfa)  +A=fa)  +A  ifa)  +  ■  +Avfa]  +  ft  is  fa).  ( A4 } 

The  ft  i  component  thus  represents  the  large-scale  trend  in 
the  data  andj^.v  may  be  the  high-frequency  noise  or  more 
isotropic  short  waves.  Figure  A I  shows  an  example  of  the 
ESD  analysis  of  the  \\  signal  of  the  DB1  data.  The  maps 
plotted  on  the  left-hand  side  from  top  to  bottom  are  f  f }  i , 
and  ftfS,  in  the  notation  of  this  appendix.  The 
scale  for  the  sequential  filtering  is  empirically  set  at  72  A\, 
36Av.  24A\,  1 2 Ay  and  6Av,  where  Av  6  m,  therefore,  the 
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Figure  A3,  (a)  Comparison  of  I  he  surface  wave  spectrum  computed  from  v-,  denoted  as  in  the 

legend,  and  the  radial  velocity,  »,  denoted  as  Su  in  the  legend,  (b)  The  ralio  between  the  wave  spectra 
outside  and  inside  the  CS. 


wavelength  of  each  11  Iter  sequence  is  432,  216,  144,  72  and 
36  nu  It  is  quite  interesting  tn  sec  that  the  cross-hatched 
patterns  of/Jj  and  J;  i  tn  the  CiS  side  of  the  current  front  arc 
enhanced  considerably  through  this  decomposition  process. 
The  choice  of  the  upper  and  lower  filter  wavelengths  (432 
and  36  in)  is  based  on  the  consideration  of  the  wave 
conditions  of  the  present  data  set  that  the  dominant  swell 
wavelength  is  about  300  m  and  the  azimuth  degrading 
becomes  serious  at  about  25  m  wavelength. 

[*:]  The  column- by -column  variance  of  the  cor¬ 
responding  signal  components  is  plotted  on  the  right-hand 
panel.  We  expect  that  the  components  and 

contain  most  of  the  swell  signal  and  the  locally- 
generated  wind  waves.  Notice  that  the  variance  of  the  wave 
signals  is  quite  small  compared  to  that  of  the  j)  \  component, 
which  is  contributed  by  a  mixture  of  large-scale  mean 
current  variations,  residual  platform  motion  and  spatial 
in  homogeneities.  The  spectra  shown  in  figure  7  are  the 
sum  of  the  individual  spectra  computed  from  the  live  wave 
signal  components  (fn^fi^  fr^/La*  and  //^l  using  v_,  In  the 
spectral  processing  of  each  decomposed  signal  component, 
a  small  area  (256.lv  x  2565a )  is  extracted  from  each  side  of 
the  velocity  front.  Detrending  of  the  area  by  high-pass  w  ith 
72 Av  (432  in)  as  the  tiller  length  parameter  is  performed 
prior  to  2D  ITT 

[u|  As  discussed  in  section  4.  the  2D  wave  spectrum  can 
also  be  calculated  from  the  sum  of  the  spectra  of  v\  and  vv 
f  igure  A 2  compares  the  wave  spectra  computed  from  these 
two  approaches  ( spectrum  of  i-r  and  the  sum  of  spectra  of  \\ 
and  \\).  The  results  arc  comparable  especially  in  the  contrast 


{ratio)  of  wave  properties  on  the  two  sides  of  the  current 
front  (Figure  A 2b). 

[ji]  The  wave  spectrum  can  also  be  calculated  from  the 
radial  velocity  directly.  I  he  transfer  function  relating  the 
radial  velocity  to  the  surface  elevation  is  [Hasselmann  and 
Hasselmann,  1 9*4 1 J 


T't  =  ^  sin  0,  j  +  i  cos  tt,J.  (AS) 

where  i\  is  the  average  incident  angle  and  k ,  the  component 
of  wavenumber  vector  in  the  radar  look  direction.  Figure  A3 
compares  the  spectra  derived  from  the  vertical  orbital 
velocity  and  the  radial  velocity.  Very  good  agreement  is 
found  in  this  swell-dominant  data  set. 


Appendix  B:  Wave-Current  Interaction 

[35]  The  wave  action  density  conservation  equation  can 
be  used  to  quantify  the  modulation  of  surface  waves  by 
surface  currents  [e.g..  Keller  and  Wright*  1 7  5 ;  Hughes, 
1978:  Phillips,  1984:  Thompson  and  Gaspamvh\  1986; 
Hwang  and  Shemdh u  1990;  Hwang,  1999.  2005), 


on_ m_ 

dt  i)t  Ox  Ot  +  Ok  Of  “ 


Uil) 


where  A  is  the  wave  action,  t  time,  a  -  (a,  y)  the  space 
vector,  k  =  (A-|,  kz)  the  wavenumber  vector,  and  Q,  source 
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terms  (the  notations  used  in  this  appendix  are  slightly 
different  from  those  used  in  the  main  text).  The  partial 
differential  equation  (Bl)  can  be  transformed  into  a  system 
of  ordinary  differentia!  equations. 


<B2) 


tix  w 


v. 


(W) 


‘‘4  -  -k[VU  -hVV.  (U4) 

tit 

where  is  the  group  velocity,  and  0  ( U<  O  the  current 

vector.  Equations  (B2)  (B3)  (B4)  can  he  solved  by  the 
method  of  characteristics  [e.g,,  Hughes,  1978;  Hwang  and 
Shemiiith  1990]  and  equation  (B4)  is  identically  the 
conservation  of  wavenumber,  equation  (II)  in  the  main 
text.  Because  ihe  wavelengths  investigated  here  are 
relatively  long,  the  magnitudes  of  the  source  functions  of 
wind  input  and  breaking  dissipation  based  on  Plant 
[1980]  or  Hughes  ||07X]  formulations,  for  example  are 
very  small.  Numerical  experiments  with  or  without  source 
functions  did  not  produce  significant  differences  in  the 
computational  results.  In  the  present  paper,  all  computations 
are  based  on  zero  external  source  terms*  The  only 
mechanism  contributing  to  the  modulation  is  the  spatial 
variation  oT  the  surface  velocity  in  the  GS,  Lor  the 
modulation  computation,  the  surface  velocity  is  assumed 
to  follow  a  1l>  bell-shaped  profile. 


where  L'i,  is  the  peak  eurrem  velocity  at  the  axis  (2  m/s 
assumed),  and  xti  a  scale  width  of  the  current  (80  kin 
assumed).  Numerical  compulations  show  that  the  level  of 
wave  variance  enhancement  attributable  to  hydrodynamic 
modulation  is  about  10  to  12  percent  lor  waves  with  initial 
(unperturbed)  wavenumbers  from  0.02  lo  0,03  rad/m.  The 
wavenumber  inside  GS  also  increases  by  a  similar 
proportion. 
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from  DOD  Ihgh  Performance  1  om puling  Modemi/ution  Program  at  the 
Naval  Research  Laboratory  Distributed  Center,  die  U.S.  Army  Engineer 
Research  and  Development  Center,  and  the  U.S.  Army  Space  am)  Missile 
Defense  C  ommand  Advanced  Research  Center,  NRI  contribution  JA 
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